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ABSTRACT 

The cross section of material in debris discs is thought to be dominated by the smallest 
grains that can still stay in bound orbits despite the repelling action of stellar radiation 
pressure. Thus the minimum (and typical) grain size Smin is expected to be close to 
the radiation pressure blowout size Sbiow Yet a recent analysis of a sample of Herschel- 
resolved debris discs showed the ratio Smin/Sbiow to systematically decrease with the 
stellar luminosity from about ten for solar-type stars to nearly unity in the discs around 
the most luminous A-type stars. Here we explore this trend in more detail, checking 
how significant it is and seeking to find possible explanations. We show that the trend is 
robust to variation of the composition and porosity of dust particles. For any assumed 
grain properties and stellar parameters, we suggest a recipe of how to estimate the 
“true” radius of a spatially unresolved debris disc, based solely on its spectral energy 
distribution. The results of our collisional simulations are qualitatively consistent with 
the trend, although additional effects may also be at work. In particular, the lack of 
grains with small Smin/Sbiow for lower luminosity stars might be caused by the grain 
surface energy constraint that should limit the size of the smallest collisional fragments. 
Also, a better agreement between the data and the collisional simulations is achieved 
when assuming debris discs of more luminous stars to have higher dynamical excitation 
than those of less luminous primaries. This would imply that protoplanetary discs of 
more massive young stars are more efficient in forming big planetesimals or planets 
that act as stirrers in the debris discs at the subsequent evolutionary stage. 
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1 INTRODUCTION 

The size distribution of dust in debris discs is set by various 
processes operating in these systems (Krivov et al. 2000). 
These include radiation pressure, collisions, transport pro¬ 
cesses, as well as mechanisms that lead to erosion of dust 
grains (Wyatt et al. 2011). Of particular importance is the 
minimum radius Smin of grains. The dust size distribution 
in most of the debris discs is steep enough for the particles 
with sizes close to Smin to be the most abundant and to carry 
most of the cross section. Therefore, they are best visible in 
scattered light, as well as in thermal emission in the near- 
and mid-infrared. In that sense, the size Smin can also be 
referred to as the typical size of grains in debris discs. 

In discs around solar- and earlier-type stars, grains 
smaller than the certain blowout limit Sbiow are expelled by 
direct radiation pressure. Thus one expects Smin to be close 
to Sbiow However, collisional simulations have shown the 
exact relation between the two to depend on the dynamical 
excitation of dust-producing planetesimals in the disc. The 
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ratio Smin/sbiow should be slightly above unity in strongly 
stirred discs (e.g., Krivov et al. 2006; Thebault & Augereau 
2007), but much larger than unity for dynamically cold 
discs (e.g., Thebault & Wu 2008; Heng & Tremaine 2010; 
Krivov et al. 2013). There are also several other parameters, 
such as the mechanical and optical properties of dust, the 
disc radius, and system’s age, which may alter Smin, Sbiow, 
and their ratio. Thus deriving the typical dust sizes from ob¬ 
servations and comparing them with the model predictions 
allows one to constrain all these disc parameters and to gain 
deeper insights into the disc physics. 

Since most of the debris discs are detected by their 
thermal emission, the easiest way to access the grain sizes 
is to analyse the temperature retrieved from the spectral 
energy distributions (SEDs). This is because the tempera¬ 
ture of different-sized grains at the same distance around 
the same star is different. Indeed, dust temperatures have 
been derived for various samples of debris discs, based on 
the IRAS, Spitzer, and Herschel data (e.g., Su et al. 2006; 
Chen et al. 2006; Rhee et al. 2007; Morales et al. 2011, 
2012; Eiroa et al. 2013; Chen et al. 2014; Kennedy & Wyatt 
2014; Mittal et al. 2015). These studies uncovered a trend of 
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temperatures increasing towards disc host stars of higher lu¬ 
minosity. Yet the interpretation of these results in terms of 
grain sizes was hampered by the degeneracy between the 
grain size and disc radius, for most of the discs in the sam¬ 
ples being unresolved. 

Booth et al. (2013) were the first to invoke a sample of 
nine resolved discs around A-type stars to break this de¬ 
generacy. They found evidence for Smin being close to Sbiow, 
as expected from the theory. Most recently, Pawellek et al. 
(2014) did a similar study for a sample of 34 resolved debris 
discs over a much broader luminosity range. They confirmed 
an increase of the minimum grain size with stellar luminos¬ 
ity, as expected from the fact the Sbiow is larger for more 
luminous stars. However, this increase was too flat to be 
consistent with Smin ~ Sbiow. Instead, the Smin/sbiow ratio 
was found to decrease from about ten for solar-type stars to 
unity for the most luminous A-type primaries — an intrigu¬ 
ing effect that needs to be explained. 

This paper presents a deeper analysis of the Smin/sbiow 
trend, shows how it can be used for unresolved discs, and 
tries to explain it. We start with an analysis of the trend, 
investigating the influence of the grain chemical composi¬ 
tion and porosity (section 2) and splitting the entire sample 
into subsamples of discs of various dustiness, radii, ages, and 
other parameters (section 3). Implications of the trend for 
the radius estimates of spatially unresolved debris discs are 
considered in section 4. Two possible explanations for the 
trend are discussed in section 5 (the surface energy con¬ 
straint on the size of the collisional fragments) and section 6 
(the possible dependence of the disc stirring level on the 
stellar luminosity). Section 7 lists our conclusions. 

2 ANALYSING THE TREND: DUST GRAIN 

PROPERTIES 

2.1 Dust compositions 

This paper uses the results of Pawellek et al. (2014) as a 
starting point. The idea was to consider a sample of well- 
resolved debris discs, for which the dust location can be 
measured from the images, and then to perform an SED 
fitting to uniquely infer the minimum grain size Smin. How¬ 
ever, Pawellek et al. (2014) assumed compact, spherical dust 
grains of pure astrosilicate (Draine 2003). The question 
arises, to what degree the results may depend on the dust 
composition. To answer it, we re-did the fitting of the entire 
sample by assuming different dust compositions. 

Since we are studying the decrease in the Smin/sbiow ra¬ 
tio with stellar luminosity, it is necessary that the blowout 
grain size exists. However, for the two M-stars of the sample 
of Pawellek et al. (2014) this is not the case for any dust 
composition. Therefore we excluded them from further in¬ 
vestigations and used only 32 objects of the previous sample. 
Another reason to exclude the M-stars was that the presence 
of the strong stellar wind of an unknown strength would in¬ 
troduce one more free parameter, making theoretical pre¬ 
dictions for Smin highly uncertain. The sample used here is 
given in Table 1. 

The SED fitting was done as follows (for details, the 
reader is referred to Pawellek et al. 2014). Many of the discs 
in the sample are suspected to have a two-component struc¬ 
ture, consisting of the main, cold outer disc (a “Kuiper belt”) 
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Figure 1. Absorption efficiency Qabs)-^) for s = 1 fim. 


and an additional, warm inner one (an “asteroid belt”). 
Since in this paper we are only concerned with the main, 
outer component, the warm one had to be looked for and, 
if present, subtracted from the SED. We did that exactly as 
described in Pawellek et al. (2014). Of the two fitting meth¬ 
ods used there, the modified blackbody emission and the size 
distribution method, here we only employed the latter one, 
since it is more reliable and gives a direct handle on the parti¬ 
cle sizes. The fitting itself was done by a simulated annealing 
algorithm, implemented in the SEDUCE code (Muller et al. 
2010). We considered five dust grain models described be¬ 
low. Altogether, we performed a second component check 
and a complete fitting of 160 objects (32 for each dust com¬ 
position), not counting additional fitting runs - e.g. to test 
extreme grain porosities. 

Specifically, we selected the compositions listed in 
Table 2 (percentages are volume fractions). We used 
the Draine (2003) optical constants for astrosilicate, the 
Zubko et al. (1996) data for carbon (their “ACAR sam¬ 
ple”) and the Li & Greenberg (1998) data for ice parti¬ 
cles. For pure astrosilicate, vacuum inclusions were added 
to simulate possible porosity. The Bruggeman mixing rule 
(Bohren & Huffman 1983) was applied to compute the re¬ 
fractive indices of mixtures, and the Mie theory was em¬ 
ployed to calculate the efficiencies. 

The dust compositions in Table 2 are ordered by the in¬ 
creasing grain temperature Tdust for a grain size of 1 /rm at a 
distance of 100 AU from a Sun-like star. For all five compo¬ 
sitions, Fig. 1 plots the absorption efficiency Qabs of 1 /rm- 
sized grains as a function of wavelength. Figure 2 presents 
the dust temperature Tdust around a star of solar luminosity 
as a function of grain radius s. The temperature curves for all 
materials are similar in shape. Small grains (s < 10 /rm) are 
warmer than blackbody, whereas the temperature of grains 
with a size larger than 10 /rm is close to the blackbody tem¬ 
perature. The grains with sizes between 0.1/rm and 1/rm 
(depending on the dust composition) are the hottest. Both 
smaller and larger grains are colder (Krivov et al. 2006). 
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Table 1. Stellar parameters sorted by stellar luminosity 


HD 

HIP 

Name 

SpT 

l/Lq 

Teft [K] 

M/Mq 

Age [Myr]-'' 

Age ref 

23484 

17439 

- 

K2V 

0.41 

5166 

0.79 

930 

1 

104860 

58876 

- 

F8 

1.16 

5930 

1.04 

200 

2 

207129 

107649 

- 

G2V 

1.25 

5912 

1.06 

2499 

1, 3, 4 

10647 

7978 

ql Eri 

F9V 

1.52 

6155 

1.12 

1307 

1, 4, 5, 6 

48682 

32480 

56 Aur 

GOV 

1.83 

6086 

1.17 

1380 

1 

50571 

32775 

HR 2562 

F5VFe+0.4 

3.17 

6490 

1.35 

449 

5, 6, 7 

170773 

90936 

HR 6948 

F5V 

3.44 

6590 

1.38 

574 

5, 7, 8 

218396 

114189 

HR 8799 

A5V 

4.81 

7380 

1.51 

71 

5, 9, 10 

109085 

61174 

rj Crv 

F2V 

5.00 

6950 

1.53 

1768 

11, 12, 13 

27290 

19893 

7 Dor 

FIV 

6.27 

7070 

1.62 

896 

10, 12 

95086 

53524 

- 

A8HI 

7.04 

7530 

1.70*' 

15 

7, 14 

195627 

101612 

01 Pav 

FOV 

7.36 

7200 

1.69 

842 

7, 8 

20320 

15197 

C Eri 

kA4hA9mA9V“ 

10.3 

7575 

1.85 

800 

12 

21997 

16449 

HR 1082 

A3IV/V 

11.2 

8325 

1.89 

44 

7, 10, 15 

110411 

61960 

(1 Vir 

AOV 

11.7 

8710 

1.91 

71 

12, 16 

142091 

77655 

K CrB 

KllVa 

12.5 

4815 

1.80" 

2345 

17, 18 

102647 

57632 

/3 Leo 

A3Va 

13.2 

8490 

1.97 

82 

12, 16, 19, 20 

125162 

69732 

A Boo 

AOp 

15.4 

8550 

2.05 

301 

12, 16 

216956 

113368 

Fomalhaut 

A4V 

15.5 

8195 

2.06 

200 

21 

17848 

13141 

u Hor 

A2V 

15.7 

8400 

2.07 

261 

6, 7, 22 

9672 

7345 

49 Get 

AIV 

16.0 

9000 

2.07 

36 

5, 7, 22, 23 

71722 

41373 

HR 3341 

AOV 

18.5 

8925 

2.16 

100 

2 

182681 

95619 

HR 7380 

BOV 

24.9 

10000 

2.33 

73 

7 

14055 

10670 

7 Tri 

AlVnn 

25.0 

9350 

2.33 

160 

12 

161868 

87108 

7 Oph 

AOV 

26.0 

9020 

2.36 

276 

7, 10, 16 

188228 

98495 

e Pav 

AOVa 

26.6 

10190 

2.37 

50 

12, 16 

10939 

8241 

q2 Eri 

AIV 

31.3 

9200 

2.47 

352 

7, 8, 24 

71155 

41307 

30 Mon 

AOV 

35.7 

9770 

2.56 

169 

12, 16 

172167 

91262 

Vega 

AOV 

51.8 

9530 

2.83 

265 

16, 25 

139006 

76267 

a CrB 

AOV 

57.7 

9220 

3.50"* 

291 

12, 16 

95418 

53910 

/3 UMa 

AlIVps 

58.2 

9130 

2.70" 

305 

12, 16 

13161 

10064 

0 Tri 

A5HI 

73.8 

8010 

4.90“* 

730 

12 


Notes: 

The effective temperatures and ages are averaged over the listed literature values. 

“Gray-Corbally notation. See App. A2 in Trilling et al. (2007) for its explanation. The majority of the stellar masses was computed 
from the luminosities by means of a standard relation M cc for main sequence stars. Exceptions for other luminosity classes or 

close binaries are the following: ^stellar mass from Moor et al. (2013) (giant); ‘^stellar mass from Bonsor et al. (2013) (subgiant); ^sum 
of the stellar masses from Kennedy et al. (2012) (close binaries); ^stellar mass from Booth et al. (2013) (subgiant). -^For each star with 
more than one age reference, the age given is the geometric mean of the values reported in those papers. 

Age references: 

[1] Eiroa et al. (2013); [2] Morales et al. (2013); [3] Lohne et al. (2012); [4] Trilling et al. (2008); [5] Moor et al. (2006); [6] Rhee et al. 

(2007); [7] Moor et al. (2015); [8] Chen et al. (2014); [9] Marois et al. (2010); [10] Chen et al. (2006); [11] Duchene et al. (2014); 

[12] Vican (2012); [13] Beichman et al. (2006); [14] Moor et al. (2013); [15] Moor et al. (2011); [16] Su et al. (2006); [17] Bonsor et al. 
(2013); [18] Bonsor et al. (2014); [19] Churcher et al. (2011); [20] Song et al. (2001); [21] Acke et al. (2012); [22] Nielsen et al. (2013); 
[23] Roberge et al. (2013); [24] Morales et al. (2011); [25] Sibthorpe et al. (2010). 


Table 2. Dust compositions 


Dust composition 

Q [g/cm3] 

50% astrosilicate -|- 50% vacuum 

1.65 

50% astrosilicate -|- 50% ice 

2.25 

100% astrosilicate 

3.30 

50% astrosilicate H- 50% carbon 

2.63 

100% carbon 

1.95 


2.2 Calculation of the blowout size 

The ratio of radiation pressure and gravitational force is 
given by (Burns et al. 1979) 

n _ ^rad _ 3-L Qpr 

^ ~ Fgrav “ IQnGMc gs ' ^ ^ 

Here, G is the gravitational constant, L the stellar lumi¬ 
nosity, M the stellar mass, c the speed of light, Qpr the 
radiation pressure efficiency averaged over the stellar spec¬ 
trum, Q the density of the grains and s the grain radius. The 
grain size that corresponds to /3 = 0.5 is the blowout grain 
size Sbiow Pawellek et al. (2014) set Qpr to one (geometric 
optics approximation). To determine Sbiow more accurately, 
we calculated Qpr for each of the material compositions as 
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Grain size [\im] 


Figure 2. Dust grain temperature Tdust(^) for ^ star with ILq 
at a radius of 100 AU. Tbb is the blackbody temperature. 


(Burns et al. 1979) 

Qpr = Qabs + Qsca(l - (C0s(l9))). (2) 

Here, Qsca is the scattering efficiency and (cos(i9)) the 
anisotropy parameter with the scattering angle '&■. 

(cos{'&)) = j f{'&) cos{'&)dQ, (3) 

where fl is the solid angle and /(i9) the phase function that 
we computed with Mie theory (Bohren & Huffman 1983). 
The averaging of Qabs and Qsca over the stellar spectra 
was done with the aid of the PHOENIX/GAIA model grid 
(Brott & Hauschildt 2005) and, for two stars with Tefi > 
10000K, ATLAS9 models (Castelli & Kurucz 2004). 

Having found Qpr(s), expression (1) should be equated 
to 0.5, yielding an equation for Sbiow. That equation has 
one root for luminous stars. For stars of moderate luminosi¬ 
ties, it may have two solutions, in which case the grains 
with sizes between the two roots are in unbound orbits. 
For the least luminous stars, the equation may not have 
solutions at all, meaning that /3 is always smaller than 0.5 
and no Sbiow exists. We solved this equation numerically. 
The results for all hve compositions and for all stars in 
our sample are shown in Fig. 3 with symbols. More exactly, 
we plot the product Sbiow(Ao/I/)'’'^'^ to eliminate the trend 
Sbiow oc L/M oc (assuming here M oc ® as appro¬ 

priate for main-sequence stars). This makes the results for 
different compositions more easily distinguishable. For com¬ 
parison, with dashed lines we depict the results calculated 
in the geometric optics approximation (e.g., with Qpr = 1) 
and assuming again M oc for all stars; multiplication 

by {Lq /renders these lines horizontal. 


2.3 Results 

Figure 4 depicts the SED fitting results for all five compo¬ 
sitions, presenting the dust temperature Tdust, the ratio of 
the true disc radius to the blackbody radius F, the grain 
size Smin, and the ratio Smin/sbiow. It demonstrates that all 
these parameters reveal more or less clear trends with the 
stellar luminosity: the dust gets warmer, the disc radius goes 
down to the blackbody value, the typical grain size increases. 



Figure 3. Blowout grain size, multiplied with [M/Mq]/{L/Lq), 
as a function of stellar luminosity. Different colours correspond to 
different dust compositions, as indicated in the legend. Symbols 
with connecting solid lines: results for size-dependent Qpr ob¬ 
tained with the Mie theory and for actual stellar masses from 
Table 1; horizontal dashed lines: results for Qpr = 1 and as¬ 
suming the main-sequence mass-luminosity relation. The verti¬ 
cal offset between the symbols and lines that represent differ¬ 
ent compositions is caused by their different bulk densities, see 
Table 2. For stars with Lq ^ T ^ IOLq , except for pure carbon, 
the second blowout limit appears in the bottom part of the plot. 
Grains in blowout orbits are those between the upper and lower 
branches. For the lowest-luminosity star in the sample (HD 23484 
with L = 0.41Lq), a single or double blowout limit only exists 
for two out of five mixtures. The strongest outliers amongst the 
symbols are stars with masses departing from the main-sequence 
mass-luminosity relation. These are the Al-subgiant /3 UMa with 
L/Lq = 58.2 and the close binaries a CrB (L/Lq = 57.7) and 
the Allll-star 0 Tri (L/Lq = 73.8). Some other apparent out¬ 
liers reflect the fact that s^iow for a-uy individual star depends on 
Qabs and Qsca averaged over the stellar photospheric spectrum 
of that particular star. For example, k CrB with LjLQ = 12.5 
is a K 1-subgiant with a temperature twice lower than that of its 
neighbouring stars in the figure. This results in the larger aver¬ 
aged values of Qabsi Qsca, and Sbiow compared to its neighbours. 
For still other stars and particular grain compositions, the stellar 
spectrum peaks at the maxima of the resonant oscillations of the 
Mie-calculated Qabs and Qsca, which also makes Sblow different 
from that of the adjacent stars. This is particularly the case for 
0 Leo with L/Lq = 13.2, if the astrosilicate-ice and astrosilicate- 
carbon mixtures are assumed. 
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Figure 4. Various dust and disc parameters as functions of stellar luminosity for different dust compositions: (a) dust temperature, (b) 
disc’s true radius to its blackbody radius, F, (c) minimum grain size Smim (d) grain size ratio Smin/sblow Similar to the previous figures, 
different colours denote different dust compositions. Symbols with error bars are fitting results for individual discs in our sample. A 
straight line of a certain colour is a best-fit trend line through the symbols of the same colour. For two objects (HD 48682 and HD 27290) 
and for two mixtures (astrosilicate-|-vacuum and astrosilicate-|-ice), the best-fit Sn,in a^nd Smin/^blow are too small for the plotting range 
of the panels c) and d). 


but the typical grain size in the blowout units decreases. All 
these results are qualitatively the same for different dust 
compositions and are consistent with those obtained previ¬ 
ously for the pure astrosilicate (Pawellek et al. 2014). 

Nevertheless, quantitative differences are apparent. In 
particular, for porous particles, the Smin/sbiow dependence 
on stellar luminosity is flatter than for compact ones. This 
is easy to understand. Figure 2 shows that for s < 3 /rm 
the temperature of porous grains is the lowest of all compo¬ 
sitions used and conversely, porous grains with sizes larger 
than 3/rm are hotter than the other mixtures. As a result. 


for early-type stars the temperature of a porous astrosilicate 
grain of size Smin is higher than that of a compact astrosili¬ 
cate particle of the same size. To reproduce the SED of a disc 
around such a star, the porous grains need to be larger than 
the compact ones. Conversely, for late-type stars a porous 
grain with size Smin is colder than a compact one and there¬ 
fore the porous particles must be smaller than the compact 
ones. The result is a steeper increase of the minimum grain 
size of porous grains with the stellar luminosity, and thus 
a gentler decrease of the Smin/sbiow ratio, compared to the 
other dust compositions. 
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factor of 



L/Lq 



Figure 5. Identification of outliers in r(L) (top panels) and Smin/sblow(f') (bottom panels) for pure astrosilicate. Left: the factors A 
and slopes B of the trend lines in the form A (L/Lq)^, calculated by excluding from the sample one data point at a time. Error bars are 
uncertainties in A and B returned by the fitting. The “cross-hairs” depict the best fit through the whole sample. The strongest outliers 
are shown in colour: red is HD 104860, magenta is rj Crv, green is Vega, and blue is /3 Tri. Right: T and size ratio versus stellar luminosity 
for pure astrosilicate with and without outliers. The colour coding is the same as on the left. The black line gives the trend line for the 
entire sample. The line of a certain colour is the best fit without the object of the same colour. 


Still, for the assumed porosity of 50%, the trend is not 
fully erased. Therefore, we tried to find the minimum degree 
of porosity that would make the Smin/sbiow ratio indepen¬ 
dent of the stellar luminosity. It turned out that for most of 
the discs the quality of the SED fits, measured by x^, gets 
poorer with increasing degree of porosity. For example, the 
disc of HD 50571 yields the reduced x^ — 1-85 for a poros¬ 
ity of zero (for pure astrosilicate). With a 50% porosity, x^ 
increases to 1.86, and with 90% to 4.55. 


We conclude that the decrease of Smin/sbiow to¬ 
wards higher stellar luminosities identified by Pawellek et al. 
(2014) for compact astrosilicate particles is pretty robust 
with respect to variation of the assumed grain properties. 


3 ANALYSING THE TREND: OUTLIERS AND 

SUBSAMPLES 

3.1 Outliers 

There is a concern about the effect of outliers in the plots 
of derived parameters versus stellar luminosity. For example, 
Tj Crv is the apparent outlier with a low disc temperature and 
F in Fig. 4a,b, and is largely the reason the best-fit lines are 
below most of the data. Also, two of the four most luminous 
stars having rather large derived Smin and small associated 
uncertainties (Vega and /3 Tri) might affect the trend lines 
of Smin and Smin/sbiow iu Fig. 4c,d. Thus some investigation 
of how much individual systems are influencing the results 
is needed to check the robustness of the conclusion of the 
Smin/Sblow trend. 

A formal way to identify the outliers and quantify their 
effect is as follows. We go over the whole sample, remove one 
data point by one, each time calculating the best-fit lines, 
and look for those discs whose removal alters the regression 
line the most strongly. We do this procedure for two plots 
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that play the major role in the rest of the paper: r(L) and 

Smin /Sbl ow (L). 

Figure 5 depicts the results for pure astrosilicate. In 
terms of F, the strongest outliers are rj Crv, HD 104860, and 
Vega. Removing t] Crv changes the slope of the trend line 
more strongly (by 0 . 12 ) than the uncertainty returned from 
the best-fit to the whole population (0.06)! For Smin/sbiow, 
three strongest outliers are P Tri, HD 104860, and Vega. This 
confirms that a few individual outliers influence the best-fit 
lines appreciably. On the other hand, dropping each of the 
objects other than those listed above has a minor effect on 
the best-fit lines. 

For these reasons, we chose to exclude three respective 
strongest outliers from the sample for the rest of the pa¬ 
per. In sections 3 and 5, where Smin/sbiow is analysed, we 
remove P Tri, HD 104860, and Vega. Besides, automati¬ 
cally excluded is the lowest-luminosity star HD 23484, for 
which — in the case of astrosilicate — no blowout limit ex¬ 
ists and thus Smin/sbiow is undefined. This implies a set of 
28 stars. In a similar style, in section 4, dealing the F-ratio 
and disc radii, we discard rj Crv, HD 104860, and Vega and 
work with a set of remaining 29 stars. 

“Refining” the sample in such a way is done in order 
not to bias the results with a few systems, which are the 
strongest outliers and for which the uncertainties are atyp- 
icially small. In doing so, we also keep in mind possible pe¬ 
culiarities of the discarded stars that may explain them be¬ 
ing outliers. For instance, HD 104860 is amongst the most 
poorly resolved sources in the sample, which makes its radius 
less certain. Besides, the low luminosity of this star makes 
us believe that the dust can be efficiently transported in¬ 
ward from the parent belt, implying that the disc radius 
measured from the image may be smaller than the true ra¬ 
dius of the parent belt. Correcting for this would return a 
smaller Smin than derived here, reducing the deviation from 
the Smin/sbiow-trend line. The strongest F-outlier, rj Crv, 
is known to be unusual in many other respects (see, e.g., 
Duchene et al. 2014). In the case of Vega, the assumed lu¬ 
minosity of 52Lq may be in error. This star is a rapid ro¬ 
tator (Peterson et al. 2006; Aufdenberg et al. 2006), which 
makes stellar parameters functions of the stellar latitude. 
The “equatorial stellar luminosity” seen by the dust could 
be as small as 28Lq and may be significantly lower than the 
“polar luminosity” of 57Lq measured from the Earth (see, 
e.g., Muller et al. 2010, for a discussion). Shifting the Vega 
point to the left in Figure 5 could make it a non-outlier. As 
far as P Tri is concerned, this is the most luminous star in 
the sample. It has some peculiarities, too, being one of only 
two giants in the sample and a close binary and, as such, 
may not be representative of the whole population. 

3.2 Extracting and comparing subsamples 

The grain sizes and the size ratios may also depend on physi¬ 
cal parameters of the systems other than the stellar luminos¬ 
ity (for instance, on disc’s fractional luminosity, disc radius, 
or system’s age). One suitable method to prove this is as fol¬ 
lows. One can split the entire sample of now 28 objects (i.e., 
without HD 23484, HD 104860, Vega, and P Tri) into a pair 
of subsamples with ni and n 2 objects (ni -|- n 2 = n = 28) 
according to the physical parameter P selected. The two 
subsamples could consist, for instance, of systems with P 



Figure 6. Size ratio vs. stellar luminosity for pure astrosili¬ 
cate, for two subsamples based on the HD numbers. Blue line 
and squares: HD number smaller than the median value 99033; 
red line and circles: HD number larger than that; black line: the 
entire sample. 

smaller and greater than the parameter’s median value Pmed, 
respectively. In that case, n\ = n 2 = 14. Then, one takes 
Smin/sbiow(E) for the two subsamples separately and derives 
the best-fit log-log trend lines with the slopes for ± SE{bi) 
and 62 ± SE{b 2 ), where SE{x) is a standard error of x. The 
final step is to find out whether the relations Smin/sbiow(D) 
in the two subsamples are statistically significantly different. 
(The null hypothesis is obviously that they are not.) To this 
end, one computes the Student’s t-score as 

t=\bi-b2\/SE{bi-b2) (4) 

with 

SE{bi - 62 ) = VSE{bi)^ + SE{b 2 )^ (5) 

and calculates the probability p that the null hypothesis is 
true on ni + n 2 — 4 = 24 degrees of freedom. We use a strict 
criterion and consider two subsamples significally different 
if the two-tailed p < 0.05. 

Since, with its 28 objects, our sample is not particu¬ 
larly large, a question arises if this procedure is feasible. 
To test this, we chose as P the HD number of the debris 
disc stars. This parameter is unphysical. Since there are no 
reasons to believe that whatever properties of debris discs 
should depend on the right assension of their primaries, our 
procedure should confirm the validity of the null hypothesis. 
The results are visualised in Fig. 6 , while the statistical pa¬ 
rameters of the regression comparison are listed in Table 3. 
The probability is p = 0.87, so that the two subsamples are 
statistically indistinguishable. 

3.3 Discs of low and high fractional luminosity 

We first take two subsamples that comprise discs with frac¬ 
tional luminosity lower and higher than the median value 
7.26 X 10~®. Figure 7 shows both the minimum grain size 
and the ratio of the minimum to the blowout size in both 
groups as a function of stellar luminosity for the reference 
material, astrosilicate. Clearly, the host stars of the discs of 
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Table 3. Comparing the regressions for Smin/^blow in pairs of 
subsamples for P < Pmed and P > Pmed 


p 

bi±SE{bi) b2±SE{b2) 

t-score 

P 

Verdict 

HD 

-0.63 ±0.09 -0.61 ±0.08 

0.17 

0.87 

same 

fd 

-0.48 ±0.11 -0.67 ±0.08 

1.40 

0.17 

same 

-^disc 

-0.53 ±0.08 -0.69 ±0.09 

1.33 

0.20 

same 

Tage 

-0.74 ±0.15 -0.58 ±0.06 

0.99 

0.33 

same 

TlOO 

-0.60 ±0.07 -0.65 ±0.10 

0.41 

0.69 

same 

Floo/Extent 

-0.70 ±0.08 -0.57 ±0.09 

1.08 

0.29 

same 

FWHM 

-0.51 ±0.15 -0.65 ±0.06 

0.87 

0.40 

same 


both classes provide a broad coverage of stellar luminosi¬ 
ties: from 1.81/0 to 581/0 and from 1.3Lq to 311/0, respec¬ 
tively. For high-/d discs, the minimum size Smin turns out 
to be Rd 5 fim, nearly independent of the stellar luminos¬ 
ity (the slope, or regression coefficient, being —0.01 ± 0.08 
only). However, Smin of low-/a discs, having the slope of 
0.17 ± 0.10), is also 2(T-consistent with being constant. The 
Smin/sbiow ratio decreases with stellar luminosity, but this 
trend is slightly stronger for high fractional luminosity discs 
(regression coefficient of —0.67±0.08) than for low fractional 
luminosity ones (—0.48 ± 0.11). Besides, the high fractional 
luminosity discs also reveal less scatter (Peason’s correlation 
coefficient of r = —O.OSIq Jg, Spearman’s Vs = —0.941q Qg) 
than those with low /a (r = —0.751 q j®, rs = —O.SOIq jg). 

Table 3 presents the formal comparison of the two sub¬ 
samples (for Smin/sbiow ratio only; but we checked that for 
Smin the results are nearly identical). The p-probability is 
0.17, so that the subsamples do not differ significantly. 

Even though the differences between the subsamples are 
of low significance, they are marginally visible. There are two 
possibilities to explain them. One is that a higher /a implies 
a more pronounced infrared excess. The SEDs of such discs 
can therefore be fitted in a more reliable way and the results 
would be less uncertain. This view is supported by the fact 
that the scatter of the data for high-/a discs around the 
trend line is smaller (i.e., the correlation between the size 
or size ratio and the luminosity is stronger). Indeed, the 
standard error of the slope for the high-/a discs (0.08) is 
smaller than that of the low-/a ones (0.11). 

Another possibility is of physical nature: it can be that 
discs with high and low fractional luminosity are physically S 

different in one or another respect. If this is true, this might 
be reflected by a systematic difference in one or another key 
parameter, such as for instance the disc radius or system’s 
age, between the two groups of discs. 


3.4 Small vs large discs and young vs old discs 

To prove whether the discs with higher and lower dust 
luminosity are physically different, we have checked the 
disc radius iodise and the age of the systems Tage in both 
groups of discs, ft turned out that the discs with /a > 
7.26 X 10“® are on the average larger (i?disc = 170 AU ± 
49 AU) and younger (Page = 488 Myr ± 685 Myr) than those 
with fd < 7.26 X 10"® (i?disc = 120AUT47AU and 
Tage = 634 Myr ± 720 Myr). Since the standard deviations of 
the radii and ages are quite large, this conclusion has to be 




Figure 7. Minimum size (top) and size ratio (bottom) vs. stellar 
luminosity, assuming pure astrosilicate, for the subsamples of low 
and high fractional luminosity discs. Blue line and squares: fd < 
7.26 X 10~®, red line and circles: fd > 7.26 X 10“®, black line: the 
entire sample. 



Figure 8. Same as Figs. 6 and 7 bottom, but for the two 
subsamples based on the disc radius. Blue line and squares: 
Tdisc < 143 AU; red line and circles: itdisc > 143 AU; black line: 
the entire sample. 


MNRAS 000, 1-16 (2015) 






































Dust grain sizes in debris discs 9 




Figure 9. Same as Figs. 6—8, but for the two subsamples based 
on system’s age. Blue line and squares: systems younger than 
284 Myr; red line and circles: systems older than 284 Myr; black 
line: the entire sample. 


Figure 10. Same as Figs. 6—9, but for the two subsamples based 
on disc’s absolute integrated brightness. Blue line and squares: 
faint discs; red line and circles: bright discs; black line: the entire 
sample. 


tested statistically. The null hypothesis is that the disc radii 
and ages in the two subsamples are indistinguishable. 

As long as the radii and ages are not distributed nor¬ 
mally, we choose the fc-sample Anderson-Darling (AD) test 
(Scholz & Stephens 1987) which is independent of a specific 
distribution and the sample size and is sufficiently sensi¬ 
tive to the tails of the distribution. We assume that the 
disc radii have an accuracy of 7% (Pawellek et al. 2014) and 
conservatively estimate the age determination accuracy as 
25% X logjQ(rage/ Myr). The AD test provides the testing 
value A^, which is given by 


m + n2 


E- 

“ rii 
1=1 


ni+n2 —1 

E 


((m +n2)Mij - juif 
j{ni +n2 - j) 


( 6 ) 


Here, Ui = 14 is the size of the i-th sample, k — 2 the 
number of samples and Mij the number of data points in 
the i-th sample smaller than the data point Zj of the pooled 


ordered sample {Zi 
standardised to 


< ... < Zn^+n^)- The A^ parameter is 


(7) 


where a is the standard deviation of A^. The null hypothesis 
is rejected at a significance level q if T > tk-i{a). For k = 2 
and a significance level oi a = 0.05, the critical percentile 
value is t\ = 1.96. For the disc radii we find T = 3.56 ± 0.79 
and for the ages T = —0.38T0.84. Thus the radii of the discs 
with high and low /u may be statistically different, whereas 
their ages are clearly not. 

Despite this result, we made two additional tests. In 
one of them we divided the whole sample into the subsam¬ 
ples of small and large disc radii, separated by the median 
radius of 143 AU (Fig. 8). In the second test, we took the 
median age of the sample (284 Myr) and divided the sample 
in young (Tage < 284 Myr) and old (Tage > 284 Myr) ob¬ 
jects (Fig. 9). The statistical parameters for the two-sample 
tests (see Eqs. 4-5) are listed in Table 3. The p-probabilities 
that the radius and age subsamples are the same are 0.20 
and 0.33, respectively. This supports the view that marginal 
differences between the discs of higher and lower fractional 


luminosity are more likely related to the difference in quality 
of the fits rather than are caused by physical reasons. 


3.5 Faint vs bright discs 

Another way to split the sample that may be illuminating 
is by the absolute integrated flux from the disc, or perhaps 
that flux divided by the disc extent. This would be a test 
for whether the results are being significantly affected by 
a systematic related to the images, since for example discs 
brighter in absolute terms may be derived to be larger sim¬ 
ply because more of the disc can be seen above the noise. 
To prove this, we now split the sample by the discs’ abso¬ 
lute integrated brightness at 100 /rm (Fig. 10) and by their 
absolute surface brightness (Fig. 11). Both faint and bright 
discs, in terms of integrated flux and surface flux, are found 
around stars of all luminosities. Judging by trend lines, we 
find no significant differences between these subsamples. The 
p-probabilities for these tests are 0.69 and 0.29, respectively 
(Table 3). 


3.6 Marginally-resolved vs well-resolved discs 

One more possibility is that trend lines are affected by the 
uncertainties with which the disc radii have been measured. 
This can be characterised by the discs’ full width at half 
maximum (FWHM) at 100 fim which can be compared to 
the Herschel/PACS point-spread function (PSF) width at 
the same wavelength, 6.8”. We split the sample into 14 discs 
with FWHM < 10.9” (i.e., FWHM < 1.6PSF) and another 
14 discs with FWHM > 10.9” (i.e., FWHM > 1.6PSF). 
The results are depicted in Fig. 12. Blue points indicating 
marginally-resolved discs are mostly found in the right por¬ 
tion of the plot and conversely, well-resolved discs concen¬ 
trate in the left part of the figure. This reveals that FWHM 
has an anti-correlation with the stellar luminosity. As dis¬ 
cussed in Pawellek et al. (2014), this reason is that more 
luminous stars in the sample are more distant on the av¬ 
erage, so that their discs often have a smaller angular size 
and are more poorly resolved. However, the trend lines in 
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Figure 11. Same as Figs. 6—10, but for the two subsamples 
based on disc’s absolute surface brightness. Blue line and squares: 
faint discs; red line and circles: bright discs; black line: the entire 
sample. 



Figure 12. Same as Figs. 6-11, but for the two subsamples based 
on disc’s FWHM at 100 /rm. Blue line and squares: marginally 
resolved discs; red line and circles: well-resolved discs; black line: 
the entire sample. 


both subsamples are pretty much the same (p = 0.40, see 
Table 3). 


4 USING THE TREND: APPLICATION TO 
UNRESOLVED DEBRIS DISCS 

We now return to the ratio of the disc radius to the black- 
body radius, F. Its dependence on the stellar luminosity was 
presented in Fig. 4b for the complete sample of 32 discs and 
different dust compositions. In view of the analysis presented 
in section 3, we now remove the three strongest outliers, 
namely rj Crv, HD 104860, and Vega. The best-fit relations 
between F and the luminosity of the central star have the 
form 

r = A {L/Lef , (8) 


Table 4. Power-law coefficients for the F ratio 


Dust composition 


A 


B 


50% astrosilicate + 50% vacuum 

5.75 

± 

0.66 

-0.40 

± 

0.04 

50% astrosilicate + 50% ice 

5.42 

± 

0.47 

-0.35 

± 

0.03 

100% astrosilicate 

8.26 

± 

1.27 

-0.55 

± 

0.04 

50% astrosilicate + 50% carbon 

6.44 

± 

0.80 

-0.41 

± 

0.04 

100% carbon 

6.61 

± 

0.74 

-0.45 

± 

0.04 


where the power-law coefficients A and B for the refined 
sample of 29 stars and all five material compositions are 
listed in Table 4. 

Table 4 shows that the trend lines of r(L) are not very 
dissimilar for all of the mixtures except for pure astrosilicate. 
In this case, the trend is the steepest and clearly predicts disc 
radii that are too large for stars of solar and subsolar lumi¬ 
nosity. This might be an argument against using pure as¬ 
trosilicate as a “reference material” for debris dust. Although 
often used, it is worth reminding that astrosilicate does not 
exist in nature. It is an artificial material that was origi¬ 
nally designed to describe the observations in the interstellar 
medium (Draine & Lee 1984), and there is no justification to 
apply it to circumstellar discs, either protoplanetary or de¬ 
bris ones. This conclusion seems also consistent with numer¬ 
ous studies of individual debris discs that showed the pure 
astrosilicate to fit the SEDs more poorly than dust composi¬ 
tions combining silicates with ice, carbon, and vacuum (e.g., 
Lebreton et al. 2012; Donaldson et al. 2013; Morales et al. 
2013; Schiippler et al. 2014, among many others). 

Relation (8) can be applied to estimate the radius of 
an unresolved debris disc. This is particularly useful, be¬ 
cause most of the known debris discs are unresolved. The 
recipe is as follows. First, one takes an SED and fits one 
or another model to it (e.g., modified blackbody emission, 
Backman & Paresce 1993) to derive the dust temperature 
Tdust. This temperature can often be found in the origi¬ 
nal papers reporting on the newly discovered excess sources, 
since such papers usually perform an SED fitting for these 
sources as well. Alternatively, Tdust can be estimated di¬ 
rectly from the wavelength Amax at which the excess emis¬ 
sion peaks: Tdust = 5100K(1 /rm/Amax). Second, one uses 
Tdust and the stellar bolometric luminosity to compute the 
blackbody radius of the disc: 


Rbb 


-idust / \-^0 / 


(9) 


Third, to get the “true” disc radius estimate, one multiplies 
Rbb by F given by Eq. (8) with coefficients A and B taken 
from Table 4 for the dust composition that is thought to be 
the most appropriate. The choice of the composition can be 
guided by the SED fitting itself. Indeed, SED fitting proce¬ 
dures usually include testing various possible compositions, 
reporting those that deliver the smallest values. 

To illustrate how the method works and to test how 
reliable it is, we can estimate the radii of a few resolved 
discs that are not in our sample, i.e., were not used to de¬ 
rive the F(T) relations in this paper. The results can then 
be comfortably compared to the actual radii, as measured 
from the resolved images. Since our sample only includes 
discs resolved by the Herschel/PACS instrument in ther¬ 
mal emission, the suitable candidates for our test can easily 
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Table 5. Resolved discs used to test the r(_L) relation. 


HD number 

L/Lq 

7dust[K] 

Rbb[AU] 

r^pred 

Rpred [ AU] 

^true[ AU] 

107146 

1.1 

51 

31 

5.5 (astrosil+vacuum) 

5.2 (astrosil+ice) 

7.8 (astrosil) 

6.2 (astrosil+carbon) 

6.3 (carbon) 

172 (astrosil+vacuum) 
163 (astrosil+ice) 

243 (astrosil) 

192 (astrosil+carbon) 
196 (carbon) 

130 

181327 

3.3 

73 

26 

3.6 (astrosil+vacuum) 
3.6 (astrosil+ice) 

4.3 (astrosil) 

3.9 (astrosil+carbon) 
3.9 (carbon) 

93 (astrosil+vacuum) 
93 (astrosil+ice) 

112 (astrosil) 

103 (astrosil+carbon) 
101 (carbon) 

89 

32297 

5.3 

83 

26 

3.0 (astrosil+vacuum) 
3.0 (astrosil+ice) 

3.3 (astrosil) 

3.2 (astrosil+carbon) 
3.1 (carbon) 

77 (astrosil+vacuum) 

79 (astrosil+ice) 

89 (astrosil) 

84 (astrosil+carbon) 

92 (carbon) 

110 


be found among the discs resolved in scattered light. We 
arbitrarily chose three such discs, namely those around a 
G2V star HD 107146 (taking the data from Ardila et al. 
2004; Williams et al. 2004; Ertel et al. 2011), an F5V star 
HD 181327 (Schneider et al. 2006; Lebreton et al. 2012), 
and an A7 star HD 32297 (Kalas 2005; Schneider et al. 2005; 
Donaldson et al. 2013). We then applied the procedure de¬ 
scribed above, separately for each of the five dust composi¬ 
tions. 


The results are listed in Table 5, and we deem them 
very reasonable. An exception is the radius estimates of the 
disc of the solar-type star HD 107146 obtained by assum¬ 
ing strongly absorbing compositions (astrosilicate, carbon, 
and their mixture). Both for HD 107146 and HD 181327, a 
better prediction for the disc radius is made by assuming a 
mixture of astrosilicate with ice or vacuum, which is consis¬ 
tent with Lebreton et al. (2012) who favour such mixtures 
over pure astrosilicate. Conversely, in the case of HD 32297, 
the discrepancy between the predicted and the observed ra¬ 
dius is the smallest for more strongly absorbing materials. A 
study over larger samples would be useful to check whether 
this is a chance effect or an indication of higher abundances 
of volatiles and/or higher porosity of dust around lower- 
luminosity stars. 


Figure 13 puts all three test objects onto the r(L) plot 
for the astrosilicate-ice mixture, showing the confidence ar¬ 
eas for the r(L) fit derived above. It demonstrates that the 
worst of our three cases, HD 32297, is a « 3(T-outlier. How¬ 
ever, black symbols show that there are a number of 3cr- 
outliers in our sample as well, so that the case of HD 32297 
is nothing extraordinary. Therefore, this example gives us a 
good grasp on the accuracy of the method. Underestimating 
or overestimating the true radius by up to a factor of two 
is possible. On any account, applying our method ensures 
a much better estimate of the true radius (89 AU-130 AU) 
than just taking the blackbody values (26 AU-31 AU for our 
three discs). 


5 EXPLAINING THE TREND: THE ROLE OF 
THE SURFACE ENERGY CONSTRAINT 


We now turn to possible explanations for the size-luminosity 
trend analysed in preceding sections. One possibility is that 
smaller collisional fragments around stars of lower luminos¬ 
ity may not be produced at all. Recently, Krijt & Kama 
(2014) pointed out that the minimum size of the collisional 
fragments created in the collisional cascade should be limited 
by the available impact energy. Denoting the fraction of ki¬ 
netic energy that is spent to create the surfaces of collisional 
fragments by rj, and the surface energy of a unit surface by 
7 , the minimum fragment radius is given by their Eq. (7): 

^ ( 10 ) 

^blow \ / / 


where 

A = (f 

~ VlOOAUy \ L J 



7 


100 erg cm“ 


( 11 ) 


and / is the mean relative velocity of colliders in the units 
of Keplerian velocity at the distance Raise from the star. 
Following Krijt & Kama (2014), / is set equal to the average 
eccentricity of the dust parent planetesimals, (e). 

Figure 14 (top) depicts the size ratio as a function of 
stellar luminosity for a subsample of 14 discs with fractional 
luminosities higher than the median value. As shown in sec¬ 
tion 3, this subsample may be more reliable than the en¬ 
tire one. The discs with high fractional luminosity reveal 
a stronger correlation of the Smin/sbiow ratio with the stel¬ 
lar luminosity, which allows an easier comparison with the 
models. 

The Smin/sbiow ratios as functions of L/Lq given by 
Eq. (10) are overplotted in Fig. 14 (top) with dashed lines 
for several values of (e). In that calculation, we assumed 
rj = 0.01 and 7 = 100 ergcm“^, with the caveat that these 
parameters are very uncertain. The regions under the dashed 
lines correspond to particles that cannot be created in col¬ 
lisions, because the impact energy would not be sufficient 
to create their surfaces. Higher eccentricities allow smaller 
grains to be produced. For (e) « 0.02-0.03, the regions ex¬ 
cluded by this constraint, together with the regions excluded 
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10 r = (4.68 + 0.69 )*(L/Lq)' 
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■0.35 + 0.06 
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Figure 13. The predicted ratio of the true disc radius to the 
blackbody radius, F, for the astrosilicate—ice mixture. The solid 
line is the best-fit relation, Eq. ( 8 ). The areas filled with dark-, 
medium-, and light-grey are Icr—, 2 ( 7 —, and 3 (t— confidence areas 
for that relation, respectively. The symbols show the actual F 
values (i.e., observed versus blackbody radius) for three selected 
test discs. Unlabeled small symbols with error bars are fitting 
results for the discs in our sample, for the same mixture. These 
are the same as blue symbols in Fig. 4b. Top: the entire sample. 
Bottom: three outliers (HD 104860, rj Crv, and Vega) excluded. 


<e>=0.01 

<e>=0.02 

<e>=0.03 

<e>=0.05 



100 


L/L(r 


7=2000 erg cm 


<e>=0.10 

<e>=0.15 



ft 


Radiation pressure blowout 


100 


L/L. 


Figure 14. The ■Smin/'^blow ratio as a function of stellar lu¬ 
minosity, assuming pure astrosilicate grains. Symbols with error 
bars: discs of our sample with /u > 7.26 X 10~®. Gray-shaded 
area: region of blowout grains. Dashed lines: original model by 
Krijt & Kama (2014), Eq. (10). Filled areas: improved model, 
Eq. (13). Grains below the dashed lines and those in filled areas 
are excluded by these models, because they should not be pro¬ 
duced. Different line and filling colours correspond to different 
degrees of the dynamical excitation, (e). Top: 7 = 100ergcm~^, 
symbols represent the full sample; bottom: 7 = 2000 erg cm“^. 


by the blowout limit, match those areas where no data points 
are present. 

As noted by Krijt & Kama (2014), there is an obvious 
possibility to further improve the model by taking into ac¬ 
count that radiation pressure should excite eccentricities of 
dust grain orbits to values higher than those of their par¬ 
ent planetesimals. Accordingly, it would be more accurate 
to attribute / to the average eccentricity of smaller solids - 
namely those, which are “immediate parents” of the colli- 
sional fragments in question. We now assume that this role 
is played by grains with radius tiSmin, where 6 > 1 is a numer¬ 
ical factor. These grains inherit the dynamical excitation of 
planetesimals, but their eccentricities are further increased 
by radiation pressure. We can adopt 


/« W(e)2 + 


l-p 


( 12 ) 
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where p is given by Eq. (1). Setting for simplicity Qpr to 
unity, which is also an assumption behind Eqs. (lO)-(ll), 
we write p = l/(2hx). Equation (10) generalises to 


/"O.Ol V 

[l 1 ^ 1 

V (e) ) 

(6)2(26*- 1)2 _ 


This is a cubic equation for x, the largest root of which is 
the Smin/sbiow ratio we are seeking. The original equation 
(10) can be obtained from (13) as a limit b —>■ oo. 

Since Eq. (13) is a cubic equation, it has two additional 
smaller roots. It is interesting to find out whether they also 
have any physical meaning. It turns out they do. The frag¬ 
ments with Smin/sbiow ratio lying between these two smaller 
roots correspond to the case where parent grains have suffi¬ 
ciently large radiation pressure-induced eccentricities to al¬ 
low production of these fragments as well. This can most 
easily be understood by looking at another limiting case, 
(e) ^ 0. In that case, Eq. (13) transforms to a quadratic 

x = 4.8x 10~^A{2bx-lf . (14) 

For the size ratios between the two roots of this equation, 
the impact energies, owing to radiation pressure-induced ec¬ 
centricities, are large enough to allow creation of fragments. 
However, for the parameter ranges considered here both 
roots of that equation are smaller than unity, i.e., corre¬ 
spond to Smin < Sblow Thus they are not physically relevant 
and are not considered further. 

The regions excluded by Eq. (13) (i.e., the ratios 
Smin/sbiow less than X, where x is the largest root of Eq. 13) 
are shown in Fig. 14 (top) as filled areas, for the same set 
of eccentricities (e). In that calculation, we set b = 10, thus 
assuming that collisional fragments are on the average ten 
times smaller than the target. (See Krivov et al. 2005, for 
justification of this choice and additional references.) 

Looking at the shape of the filled regions shown in 
Fig. 14 (top), the refined model may seem to reproduce the 
trend in the data points more poorly than the original model. 
However, a closer similarity between the model constraints 
and the data can be achieved by varying the coefficient A in 
Eq. (13). This is demonstrated by Fig. 14 (bottom) where 
we increased the coefficient A (see Eq. 11) by a factor of 20. 
The latter may correspond, for instance, to a 20 times higher 
surface energy per unit surface area, 7 (still not unrealistic, 
e.g., for icy material) or to a 20 times lower energy fraction 
that goes to the surface creation, 77 (not unrealistic either). 
With this choice, a reasonable match to the trend seen in the 
data is achieved for somewhat higher values of (e). Indeed, 
one sees that for (e) « 0.10, the region in Fig. 14 (bottom) 
excluded by our model could masquerade as a decrease in 
Smin/sbiow Seen in the data. 


6 EXPLAINING THE TREND: THE ROLE OF 
THE STIRRING LEVEL 

In this section, we check another possible explanation for the 
trend in grain sizes with the stellar luminosity. The mecha¬ 
nism that we will address here is dissimilar, albeit not com¬ 
pletely unrelated, to the microphysical one discussed in the 
previous section. It is associated with the balance between 
the production and loss rate of small grains, controlled by 
the stirring level of larger bodies. 


6.1 Idea 

Thebault & Wu (2008) predicted the grain size distribution 
to depend on the degree of dynamical excitation of the dust- 
producing planetesimals, (e). If the planetesimals have a low 
dynamical excitation which, however, is still high enough 
for collisions to be mostly destructive, then the low colli¬ 
sion velocities between large grains, that are not suscepti¬ 
ble to radiation pressure, would decrease the rate at which 
small grains are produced. However, the destruction rate of 
these small grains is set by eccentricities induced by radia¬ 
tion pressure and remains the same. This should result in a 
dearth of small dust. The maximum of the size distribution, 
i.e. the ratio Smin/sbiow, would shift to larger values. This 
is easy to quantify. Following the explanation above, the 
size distribution should peak at grain sizes Smin, for which 
the radiation pressure-induced eccentricity, /?/(! —/3), equals 
the eccentricity inherited from the planetesimals, (e). Since 
p « O.5(sbiow/Smin) (assuming Qpr = 1), this gives 

Smin/Sblow « ((e)”^ -k l) /2. (15) 

However, Eq. (15) is just a rough estimate. Obviously, 
the effect has to be confirmed by collisional simulations that 
include a more realistic treatment of collisional and radiation 
pressure forces. Such a study is still missing, since the paper 
by Thebault & Wu (2008) was confined to discs of A-stars, 
and the size distribution was computed with a collisional 
code that did not include cratering collisions, rebounds and 
sticking. 

6.2 ACE runs 

To quantify the effect more accurately, we performed four 
runs of our ACE code (Krivov et al. 2013), probing two cen¬ 
tral stars (A2V with 17.4Lq and G2V with ILq) and two 
typical stirring levels (average planetesimals’ eccentricity of 
0.1 and 0.01, average inclination according to the energy 
equipartition). The simulations included stellar gravity, di¬ 
rect radiation pressure and Poynting-Robertson drag, and a 
wealth of possible collisional outcomes (disruptive, cratering, 
rebounding, and sticking collisions). In all of the cases, the 
initial disc mass was taken to be 30M@ (in the bodies of up 
to 100km radius), in order to arrive at the typical dust frac¬ 
tional luminosity level in our sample. The disc radius was set 
to 100 ± 10 AU. We also made a number of other standard 
assumptions. In particular, we assumed compact astrosili- 
cate from Draine (2003) as a material composition, the crit¬ 
ical fragmentation energy from Benz & Asphaug (1999), etc. 
Each disc was evolved until a quasi-steady state, as defined 
in Lohne et al. (2008), was reached. 

6.3 Results 

The resulting size distributions are shown in Fig. 15. They 
confirm the effect of the maximum in the size distribution 
shifting towards larger values for discs with lower dynam¬ 
ical excitation. However, the size distributions simulated 
with ACE cannot be closely approximated by power laws 
with a sharp lowest cutoff. In the case of G2-star discs 
(sbiow = 0.46 /rm), the position of the maximum in the size 
distribution, 2.6 /rm for (e) = 0.1 or 10 fim for (e) = 0.01, can 
still be taken as Smin- In the discs around an A2-star, the size 
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Figure 15. Simulated size distributions in the fiducial discs 
around A2-stars (thick lines) and G2-stars (thin lines) with higher 
(red) and lower (blue) level of stirring. Plotted is the size distribu¬ 
tion in the parent ring, assumed to be located at 100 AU. Vertical 
dashed lines mark the radiation pressure blowout limit of A2-stars 
(thick) and G2-stars (thin). Note that two blowout values exist 
for the G2-stars. The grains in unbound obrits are those between 
these two. 



Figure 16. Symbols with error bars represent the discs of our 
sample (assuming pure astrosilicate) and the thick black line is the 
best-fit through these data. Asterisks mark approximate values of 
^min/sblow extracted from ACE simulations, as described in the 
text. The pairs of asterisks corresponding to the same dynamical 
excitation, (e), are connected with dashed lines for illustrative 
purposes. 


distribution develops a “plateau” between Sbiow = 3.0 fim 
and a shallow maximum at 12 /im (for (e) = 0.1) or 50 ^m 
(for (e) = 0.01). The different-sized grains in these ranges 
contribute to the cross section almost equally. In that case, 
we take a geometric mean between the two ends of the 
plateau as a proxy for Smin. 

Figure 15 also shows that for the G2V stars, two 
blowout limits exist. The “sub-blowout grains” of sizes 
smaller than <0.1^m are also is bound orbits and make 


a large contribution to the total cross section of dust. How¬ 
ever, their absorption efficiency is much lower than that of 
the large bound grains. Accordingly, their contribution to 
the observed emission is minor, at ~ 10% level at all wave¬ 
lengths. For this reason, the sub-blowout grains are not di¬ 
rectly relevant to the above discussion of the lower cutoff in 
the size distribution. Nevertheless, these grains still implic¬ 
itly affect the size distribution of large bound grains, since 
collisions with them “erode” the population of particles with 
sizes just above the main blowout size. 

Figure 16 compares the resulting Smin/sbiow ratios found 
with the ACE runs with our sample. As expected, lower (e) 
lead to larger Smin/sbiow, and Smin are roughly consistent 
with Eq. (15). What was not really expected prior to our 
modelling, however, is that the lines of constant dynamical 
excitation are tilted. The collisional modelling shows that for 
stars of earlier spectral types, Smin becomes closer to Sbiow 
even if (e) = const — an effect that was not previously re¬ 
ported in the literature. Coincidentally or not, this is exactly 
the qualitative trend we found in the observational data. 

However, the observed trend is stronger, and is not con¬ 
sistent with a constant (e) assumption for stars across the 
full luminosity range. Instead, the data seem to be best re¬ 
produced by assuming (e) between 0.01 and 0.1 for G-stars 
and (e)>0.1 for A-stars. This impression is supported by 
the detailed collisional modelling of several individual ob¬ 
jects done previously. For instance, Lohne et al. (2012) and 
Schiippler et al. (2014) find (e) to be at the level of a few per 
cent for the discs of a GOV-star HD 207129 and a K2V-star 
HIP 17439, respectively, whereas higher values (e) ~ 0.1 are 
favoured for the discs of an A5V-star HR8799 (where it is 
expected based on the presence of a substantial outer halo 
of small grains, see Matthews et al. 2014) and an AOV-star 
Vega (Muller et al. 2010). We stress, however, that these 
claims should not be overinterpreted. We are only discussing 
general statistical trends. Individual systems show an appre¬ 
ciable scatter, and there are certainly discs that do not follow 
those general trends. 

There is one more aspect related to the stirring level. 
If the discs around more luminous stars are excited more 
strongly, they should possess more pronounced halos. In that 
case, the disc radii of early-type stars retrieved from the im¬ 
ages may be overestimated, and so, the dust grain sizes may 
be underestimated (Pawellek et al. 2014). Gorrecting for this 
effect would push the data points on the right of Fig. 16 
up, flattening the Smin/sbiow dependence on stellar luminos¬ 
ity. However, this would make the excitation level of these 
discs estimated from the plot lower than before the correc¬ 
tion. This would imply the weaker halos, pushing the data 
points of luminous stars back down and resulting in smaller 
Smin/sbiow at high luminosities and a steeper dependence. 
Thus, the Smin/sbiow(T) dependence plotted in Fig. 16 may 
be, so to say, “self-regulating” and should be relatively ro¬ 
bust with respect to the uncertainties in measuring the disc 
radius from the images. 

Another remark is that the strirring models considered 
here and the surface energy models addressed in section 5 
are not independent. If the stirring level does depend on 
stellar luminosity as proposed here, this will affect the sur¬ 
face energy constraint discussed in section 5, changing the 
filled areas shown in Fig. 14. The best way to address this 
would be to combine both effects in a single model. This 
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can be done in the future by implementing the surface en¬ 
ergy constraint directly into the kernel of the collisional code 
that controls the production of fragments (the ACE simu¬ 
lations presented here did not include the surface energy 
constraint). 

Can we expect the discs of more luminous stars 
to be more strongly stirred on the average than those 
of low-luminosity stars? In principle, yes. Indeed, it is 
known that the submillimeter dust masses in protoplane¬ 
tary discs, the progenitors of debris discs, are roughly pro¬ 
portional to the masses of their central stars (see, e.g.. 
Fig. 5 in Williams & Cieza 2011, and discussion therein). 
It is possible that discs of more massive (or more lu¬ 
minous) stars, being more massive, were more “success¬ 
ful” in building large planetesimals that may later act 
as stirrers for debris discs (the so-called “self-stirring sce¬ 
nario”, e.g., Wyatt 2008; Kennedy & Wyatt 2010). This is 
directly supported by planetesimal accretion simulations 
(e.g., Kenyon & Bromley 2008). Alternatively or addition¬ 
ally, such protoplanetary discs may have formed more giant 
planets, and/or these planets may have experienced more 
vigorous migration or scattering in the past, which might 
have also resulted in a higher dynamical excitation of the 
debris discs emerged in these systems (the “planetary stir¬ 
ring scenario”, e.g., Mustill & Wyatt 2009). Indeed, obser¬ 
vational (e.g., Johnson et al. 2010; Reffert et al. 2015) and 
theoretical work (e.g., Ida & Lin 2005; Kennedy & Kenyon 
2008; Alibert et al. 2011; Mordasini et al. 2012) find that 
giant planets are more frequent around more massive stars, 
and that those planets are typically more massive. Further 
work is required, however, to validate or falsify these possi¬ 
bilities. 


7 CONCLUSIONS AND DISCUSSION 

Analysing a sample of Herschel-resolved debris discs, 
Pawellek et al. (2014) found the ratio of the minimum (or 
typical) grain size Smin to the radiation pressure blowout 
size Sbiow to systematically decrease with the increasing lu¬ 
minosity of the central stars. Here we investigate how robust 
the trend is and attempt to find possible explanations for it. 
Our conclusions are as follows: 

(i) The decrease of Smin/sbiow with increasing luminos¬ 
ity of the central stars persists, no matter which material 
compositions and porosity of dust grains are assumed. The 
trend is gentler for porous grains, but can be completely 
erased only for unrealistically high porosities. 

(ii) The minimum grain size itself is consistent with being 
constant, Smin ~ 5 pm ± 0.3 pm, across the full luminosity 
range of the sample. 

(iii) We have tested the subsamples of discs with lower 
and higher fractional luminosity, smaller and larger radii, 
younger and older ages, lower and higher absolute integrated 
flux, lower and higher absolute surface brightness, as well 
as those that are marginally and well-resolved. In terms of 
Smin/sbiow(L), all of these subsamples are found to be statis¬ 
tically indistinguishable. Marginal differences are only visi¬ 
ble between the discs with higher and lower fractional lumi¬ 
nosity in our sample. We argue that these are likely related 
to SED-fitting aspects rather than reflect any real physical 
difference between the dustier and the less dusty discs. 


(iv) As a by-product of this study, we derive empirical 
formulae for the ratio T of the true disc radius to its black- 
body radius, as a function of stellar luminosity. Since the 
blackbody radius can easily be found from the SED of any 
excess source, this offers a recipe of how to estimate the true 
radius of unresolved debris discs, thus breaking the notori¬ 
ously known degeneracy between the dust grain sizes and 
the dust location. 

(v) The results of our collisional simulations reproduce 
the Smin/sbiow trend with the stellar luminosity seen in the 
sample. 

(vi) Additional effects may also be at work, contribut¬ 
ing to the observed trend. For instance, the surface energy 
constraint on the size of the smallest collisional fragments 
identified by Krijt & Kama (2014) may be responsible for 
the absence of grains with small Smin/sbiow ratios in discs of 
lower luminosity stars. 

(vii) It is also possible that the Smin/sbiow-trend is related 
to the degree of stirring of the dust-producing planetesimals. 
Indeed, a better agreement between the data and the col¬ 
lisional simulations can be achieved by assuming that the 
discs of earlier-type stars are more strongly excited (on the 
average) than those of later-type stars. This would imply 
that protoplanetary discs of more massive young stars are 
more efficient in forming big planetesimals or planets that 
act as stirrers in the debris discs at the subsequent evolu¬ 
tionary stage. 

Regardless of whether the observed trend is related to 
microphysics of collisional dust production or stirring levels 
of planetesimals, there is a close relation between these two 
types of models. Both involve small grain dynamics being set 
by radiation pressure. This emphasizes again that dust grain 
dynamics is a key ingredient in debris disc models, without 
which confronting simulations and observations would be 
unthinkable. 

In principle, it is possible to find alternative explana¬ 
tions for the trend by allowing various parameters of dust 
and/or planetesimals to systematically vary with the type of 
the central stars. These may include, for instance, the chemi¬ 
cal composition, degree of porosity, or tensile strength. How¬ 
ever, such an assumption is not substantiated by any pieces 
of observational evidence. Nor is it supported by theoretical 
expectations. For example, there is no compelling evidence 
for a systematic difference in the chemical composition in 
the cold parts of the discs of Herbig Ae/Be stars and those 
of T Tan stars, which are progenitors of debris discs around 
early and later-type stars, respectively (e.g., Dutrey et al. 
2014). 
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